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二 层 分 解 技 术 在 电价 预测 中 的 应 用 研究 
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(广东 工业 大 学 , 广州 510006) 


摘 要 : 针对 电价 波动 幅度 大 且 预 测 精度 低 的 问题 ， 提 出 了 二 层 分 解 技 术 与 神经 网 络 组 成 的 电价 多 步 预 测 模型 。 该 模 
型 首先 采用 集合 经 验 模 式 分 解 将 原始 电价 序列 分 解 为 一 系列 分 量 ， 变 分 模 态 分 解 将 第 一 层 分 解 产 生 的 最 高 频率 分 量 进 
一 步 分 解 为 一 系列 模 态 分 量 ,所 有 分 量 采用 神经 网 络 模型 进行 预测 ,并 使 用 纵横 交叉 算法 对 神经 网 络 的 参数 进行 优化 ， 
最 后 登 加 所 有 子 序 列 ， 得 出 预测 电价 值 。 仿 真 结果 表明 ， 所 提出 的 模型 相 比 其 他 混合 模型 具有 更 好 的 预测 性 能 ， 且 
用 价值 高 。 
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Research on two-layer decomposition technique for predicting electricity price 


Yin Hao, Zeng Yun, Huang Shengquan, Dong Zhen 
(College of Automation Guangdong University of Technology, Guangzhou 510006, China) 


Abstract: To solve the problems that the electricity price has a lot of volatility and a low accuracy , this paper proposesa multi- 
step prediction model composed of two-layer decomposition technique and back propagation neural network(BP) . Firstly, 
variational mode decomposition(VMD) is specifically applied to further decompose the high frequency intrinsic mode functions 
(IMFs) generated by ensemble empirical mode decomposition(EEMD) into a number of modes. Next, BP model optimized by 
crisscross optimization algorithm(CSO) algorithm is utinized to forecast all sub-sequences, and the forecast series of electricity 


price is obtained by adding up all sub-sequences. The simulation result shows that the proposed model has superior performances 


for other hybrid models , and it has a great practical values. 


Key Words: two-layer decomposition; crisscross optimization algorithm; multi-step prediction, neural network; electricity price 


prediction 
0 引言 电价 受 很 多 复杂 因素 影响 ,因此 考虑 到 用 因果 关系 预测 法 具有 
一 定 的 难度 ,所 以 许多 研究 人 员 转 向 使 用 时 间 序 列 法 进行 电价 
能 是 当今 世界 必 不 可 少 的 绿色 无 污染 能 源 , 在 人 们 日 常 预测 四。 最 常用 的 时 间 序 列 法 分 为 : 统计 模型 器 、 混 合 模型 中、 
生活 中 扮演 着 重要 角色 。 在 电力 市 场 竞 争 越 来 越 激烈 的 情况 下 ， ”人 工 智能 模型 。 而 现 已 有 大 量 人 工 智能 法 应 用 于 不 同 的 预测 领 
电力 可 以 像 货物 一 样 在 市 场 环 境 下 自由 交易 ,， 因 此 电价 在 很 大 域 ， 如 人工 神经 网 络 中 、BP、 支 持 向 量 机 加 〈support vector 
程度 上 反映 了 电力 市 场 的 供求 关系 。 由 于 电价 受 负荷 、 天 和气、 machine，SVM) ， 如 文献 [9] 使 用 粒子 群 算法 (particle swarm 
商务 交易 、 日 常 活动 以 及 发 电 侧 报价 等 人 为 因素 的 影响 ， 使 得 optimization，PSO) 对 神经 网 络 的 闵 值 和 权 值 进行 优化 ， 来 提 


EE 价 具有 随机 性 、 不 稳定 性 、 非 线性 的 复杂 特点 中， 故 建立 精 高 神经 网 络 的 预测 能 力 ， 并 提出 了 PSO-BP 的 边际 电价 预测 模 
准 的 电价 预测 模型 面临 着 巨大 的 挑战 。 电价 波动 大 小 与 电力 市 型 。 将 优化 的 神经 网 络 模型 和 未 优化 神经 网 络 模型 的 计算 结果 
场 风 险 紧 密 相关 争 ， 因 此 ， 准 确 的 电价 预测 对 于 电力 系统 和 电 进行 比较 分 析 ， 可 知 优化 后 的 模型 预测 结果 和 实际 数据 相差 更 
力 市 场 具 有 重要 的 意义 。 、。 对 神经 网 络 参数 优化 是 一 个 大 规模 多 峰 优 化 问题 ,PSO 算 
电价 预测 方法 分 为 二 类 ,一 类 是 通过 寻找 自 变量 与 预测 值 法 在 优化 时 存在 收敛 速度 慢 、 易 陷入 局 部 最 优等 丈 端 。 
因果 关系 的 因果 关系 预测 法 钻 ， 男 一 类 是 时 间 序 列 法 ， 此 方法 由 于 很 多 数据 系列 具有 非 平稳 性 和 非 线性 的 特点 ， 如 电价 、 
是 假设 预测 值 和 历史 电价 相关 联 。 如 上 所 述 ， 在 实际 生活 中 ， 风速 、 负 荷 等 ， 然 而 单一 预测 模型 已 不 能 精准 预测 复杂 的 数据 
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系列 。 因 此 很 多 学 者 一 直 努 力 使 用 不 同 的 分 解 技术 去 处 理 数据 absolute percentage error，MAPE) 判断 最 佳 分 解 个 数 ， 如 图 3 
系列 中 存在 的 复杂 问题 ， 以 便 更 精准 的 进行 预测 。 如 文献 [10] 所 示 , 是 MAPE 随 着 分 解 个 数 变化 的 曲线 , 随 着 分 解 个 数 增 加 ， 
提出 一 种 小 波 包 变 换 (wavelet packettransform，WPT) 和 非 参 MAPE 逐渐 减 小 ， 但 是 在 分 解 个 数 大 于 8 的 时 候 ， 预 测 误差 稍 
数 (generalized auto regressive conditional heteroskedasticity ， 微 有 点 增加 ， 所 以 本 文 模型 中 选择 EEMD 分 解 个 数 为 8， 分 别 
GARCH) 混合 模型 对 日 前 电价 进行 预测 ， 实 验 表明 ，WPT 能 是 IMF1、IMF2、…、IMF7 和 一 个 剩余 量 。 如 图 4 所 示 ,IMFs 
很 大 程度 提高 预测 精度 。 文 献 [11] 提 出 一 种 基于 小 波 变 换 频率 从 高 到 低 排 列 , 每 个 IMF 反映 了 存在 于 序列 中 的 不 同 振荡 
(wavelet transform, WT) 和 遗传 算法 (genetic algorithm, GA) 模式 ，IMF1 是 反映 电价 详细 信息 的 最 高 频率 信号 ， 剩 余 量 是 
天 化 SVM 新 型 混合 模型 。 尽 管 现在 很 多 单 分 解 技术 已 广泛 运 频率 最 低 的 信号 ， 显 示 了 电价 波动 的 基本 趋势 。 

用 于 预测 问题 上 ， 如 文献 [12,13] 中 分 别 使 用 了 经 验 模式 分 解 
(empirical mode decomposition, EMD ) 和 集合 经 验 模式 分 解 。 me | 


一 一 昆士兰 州 


然而 高 频率 本 征 模 函数 (intrinsic mode functions ，IMFs)， 特 别 2 
是 IMF1， 很 容易 受到 许多 随机 因素 的 影响 鸣 ， 从 而 加 大 了 模 200 
拟 难 度 。 此 外 ， 本 文 发 现 由 于 很 多 非 线 性 和 非 稳 定性 的 原始 序 Si% 
列 导 致 了 更 多 的 高 频率 IMFs 产生 55， 所 以 一 层 分 解 技术 不 能 100 
td 

为 了 解决 上 述 问 题 ， 结 合 单 分 解 技术 的 优势 ， 本 文采 用 了 "0 400 600 700 
二 层 分 解 技术 ， 提 出 ] 站 于 集合 经 凤 模 式 分 角 变 分 模 态 分 解 9 

到 1 二 个 电力 市 场 的 原始 电价 序列 

和 纵横 交叉 算法 优化 的 神经 网 络 组 成 的 新 型 混合 预测 模型 
(EEMD-VMD-CSO-BP ) 来 实现 电价 的 多 步 预测 。EEMD 将 原 
始 电价 序列 分 解 为 一 系列 频率 由 高 到 低 的 IMFs 和 一 个 剩余 量 
CResidual ) ， 为 了 减少 数据 序列 的 非 稳 定性 ，VMD 进一步 将 
高 频率 IMFs 分 解 为 多 个 模 态 ， 并 用 纵横 交叉 算法 (crisscross 
optimization algorithm，CSO ) 优化 神经 网 络 的 权 值 和 浆 值 来 提 
高 预测 精度 。 本 文 分 别 以 以 澳大利亚 的 新 南 威尔士 州 和 昆 士 南 
州 的 电力 市 场 某 月 的 电价 数据 作为 样本 ， 进 行 了 一 步 、 三 步 、 
五 步 电价 预测 ， 并 对 实验 结果 进行 了 分 析 。 根 据 对 比 模型 仿 
结果 表明 ， 本 文 提出 的 混合 模型 预测 精度 高 于 其 他 混合 模型。 oon |. lee, 
1 EEMD-VMD-CSO-BP 混合 模型 a 

本 文 仿真 数据 分 别 来 源 于 澳大利亚 新 南 威尔士 州 ( 电 力 市 EE 
场 1) 和 昆 士 南 州 《电力 市 场 2) 的 电力 市 场 实测 电价 数据 ， 由 
于 不 同 区域 有 着 不 同 的 规模 ， 如 人 口 、 地 理 位 置 、 产 业 结 构 和 项 测 电价 序列 
气候 特征 ,这些 存在 的 不 同 因素 导致 了 电价 波动 出 现 了 很 大 的 

图 2 EEMD-VMD-CSO-BP 混合 模型 流程 图 

差异 ， 因 此 本 文采 用 二 个 电力 市 场 的 电价 可 以 全 面 系统 评估 所 
提 模 型 的 有 效 性 和 实用 性 。 本 文 以 每 半 小 时 为 1 个 数据 点 ， 每 b) 显 而 易 见 ,在 图 5 中 ,除了 具有 高 频率 振荡 特性 的 IMF1， 
天 共有 48 个 观测 值 ， 根 据 这 二 个 电力 市 场 某 月 上 旬 记 录 电 价 含有 很 多 噪声 序列 ， 其 余 IMFs 和 剩余 量 波动 较 平稳 ， 且 易 能 
值 (共有 720 个 观测 值 ) ， 可 以 得 到 图 1 所 示 新 南 威尔士 州 和 准确 的 预测 ， 由 于 IMF1 会 增加 预测 难度 ， 为 了 使 预测 效果 更 
昆 士 南 州 的 电价 序列 ， 并 利用 原始 序列 前 600 个 数据 作为 训练 好 ， 在 第 二 层 分 解 中 ， 利 用 VMD 将 IMF1 分 解 为 多 个 模 态 分 
样本 ， 对 第 600~700 个 电价 数据 进行 预测 。 量 ， 根 据 本 文 所 提 模 型 来 大 致 判断 VMD 分 解 个 数 ， 如 图 4 所 

本 文 所 提出 EEMD-VMD-CSO-BP 混合 模型 基本 框架 如 图 示 ， 当 分 解 个 数 趋 于 9 后 MAPE 处 于 相对 平稳 状态 , 所 以 大 致 


2 所 示 ， 所 提 混 合 模型 主要 框架 是 基于 二 层 分 解 技术 机 制 ， 以 选择 分 解 个 数 为 9, 分 别 是 model、mode2、…、mode8、mode9， 

下 提 及 实验 仿真 都 是 以 电力 市 场 1 的 半 小 时 数据 为 例 。 详 细 步 如 图 6 所 示 。 

又 如 下 ; c) 采 用 纵横 交叉 算法 优化 的 神经 网 络 模型 对 所 有 子 序列 进 
a) 在 第 一 层 分 解 中 ， 为 了 降低 电价 序列 非 线 性 和 非 平稳 性 行 预测 ， 子 序列 包括 model、mode2、…、mode9、IMF2、……、 

的 特点 ，EEMD 用 于 将 原始 电价 序列 分 解 为 多 个 分 量 ， 可 大 概 IMF7、 剩 余 量 。 闪 加 全 部 子 序列 预测 值 ， 所 得 值 即 实际 预测 电 

根据 EEMD-CSO-BP 混合 模型 的 平均 绝对 百分比 误差 (mean 价值 。 
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EEMD 分 解 层 数 
图 3 预测 误差 与 EEMD 分 解 层 数 变 化 图 (电力 市 场 1) 
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图 4 ”预测 误差 与 VMD 分 解 层 数 变 化 图 (电力 市 场 1) 


1 1 
7 8 9 10 


Mode9 Mode8 Mode7 Mode6 Mode5 Mode4 Mode3 Mode2 Model 


0 100 200 


300 400 
Time(halfhour) 


500 600 700 


5 ”EEMD 分 解 电力 市 场 1 的 原始 电价 数据 
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2 ”二 层 分 解 技 术 


2.1 集合 经 验 模式 分 解 

经 验 模 式 分 解 (empirical model decomposition, EMD ) 是 一 
种 自 适 应 信号 时 频 处 理 方法 ， 适 用 于 将 非 线性 、 非 平稳 的 复杂 
信和 号 分 解 为 一 些 频率 由 高 到 低 的 本 征 模 函数 和 一 个 剩余 量 。 集 
合 经 验 模式 分 解 (EEMD) 是 在 经 验 模 式 分 解 (EMD) 的 基础 上 
进行 改进 ， 通 过 加 入 高 斯 白 噪声 这 个 变量 ,来 处 理 经 验 模 式 分 
解 中 常 遇 到 的 模 态 混 辣 现象 ， 并 最 大 化 保存 了 原始 信号 。 本 文 
EEMD 分 解 中 里 的 三 个 重要 参数 : 白 噪 声 久 、 整 合 数 n、 复 制 
次 数 M 分 别 设置 为 0.4、9、200。EEMD 的 详细 原理 参考 文献 
[16]， 具 体 分 解 步骤 如 下 : 

a) 在 原始 数据 序列 中 加 入 高 斯 白 噪声 。 

b) 将 含有 白 噪声 的 EMD 序列 分 解 为 多 个 IMF 和 一 个 剩余 


c) 通 过 加 入 不 同 的 白 噪声 来 重复 步 又 (1) 和 (2), 并 得 到 相应 
的 IMF 和 剩余 量 ， 重 复 次 数 即 整合 数 n。 

gd 将 所 有 的 IMF 的 均值 和 所 有 剩余 量 部 分 均值 作为 最 终 
结果 。 
2.2 变 分 模 态 分 解 

变 分 模 态 分 解 (VMD) 是 DRAGOMIRETSKIY 在 2014 年 
时 出 的 一 种 新 的 非 递归 信号 处 理 技术 07, 能 够 自 适应 将 实 值 信 
号 分 解 成 一 些 拥有 特定 稀 玻 属性 的 离散 模 态 。VMD 实质 上 是 
一 个 变 分 问题 ， 为 了 处 理 该 问题 ， 可 以 将 所 有 模 态 看 做 拥有 不 
用 中 心 频 率 的 有 限 宽带 ， 利 用 交 蔡 方向 乘 子 法 ， 逐 步 将 所 有 模 
态 调 解 到 对 应 的 基 频 带 ， 最 后 得 到 各 模 态 及 各 模 态 对 应 的 中 心 
频率 。 为 了 评估 每 个 模式 的 频 宽 ， 应 该 考虑 以 下 步 又: 

a) 对 于 原始 信号 f(?), 利用 Hilbert 变换 进行 时 频 处 理 , 进 
而 获得 每 个 模 态 函数 w(t) 的 解析 信号 ， 以 及 u(t) 的 单 边 频谱 。 
b) 为 了 使 每 个 模 态 的 频谱 调制 到 基带 上 , 则 在 w(t) 中 混合 
了 相应 的 指数 项 ， 使 得 调整 到 所 估计 的 中 心 频率 。 

c) 使 用 以 上 已 调解 信号 的 高 斯 平滑 度 H1 来 计算 每 个 模 态 
的 带宽 。 因 此 ， 得 到 一 个 变 分 约束 问题 : 

| (1) 


{ux 记 by ol (oO+ 去 e 0 
好 是 分 解 后 获得 的 KK 


st = f(t) 
个 模 态 分 量 ,其 中 好 = {tux} 他 } 为 每 个 模 态 对 应 的 中 心 


— 


在 式 (1) 中 ，f() 是 原始 数据 信号 


频率 ,其 中 fs = {wi wk k=1,2,3…K 。 考 虑 到 惩罚 项 


和 拉 格 朗 日 乘 数 K， 可 以 将 上 式 转 换 为 无 约束 问题 ， 形 成 了 扩 
展 拉 格 朗 日 表达 式 ， 如 式 (2) 所 示 。 
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其 中 : a 为 数据 保 真 度 约束 的 平衡 参数 ， 可 以 在 噪声 存在 的 情 
况 下 也 能 保持 重 构 的 准确 性 ， X(t) 是 拉 格 朗 日 算 子 。 通 过 
ADMM 方法 在 每 次 迭代 优化 中 可 以 找到 上 式 表 达 式 的 鞍点 [3]， 
也 就 是 式 (D) 中 的 最 优 解 。 其 中 帮 和 wi 的 表达 式 如 下 : 


、 A(w 
-ar 


sntl 放大 
人 ty 4 2a(w— w 了 

[wha (wj am 
wr 0 


[hi& (wj aw 

0 

其 中 : w 为 信号 频率 , n 表示 办 代 次 数 ，f(w)、 训 、 部 "(w)、 
Aw) 分 别 是 了 QO) 、w()、tw*) 、X(E) 的 傅 里 叶 变 换 。 

3 ”纵横 交叉 优化 算法 


纵横 交叉 算法 (CSO) 是 由 横向 交叉 和 纵向 交叉 两 种 算 子 
组 成 。 每 次 迭代 过 程 中 ， 都 会 进行 横向 交叉 和 纵向 交叉 ,交叉 


3.2 ”纵向 交叉 操作 


纵向 交叉 操作 为 了 使 陷入 局 部 最 优 的 粒子 跳出 来 ， 而 又 


干扰 其 他 维 的 粒子 。 


于 不 同 维 所 处 范围 不 同 ， 因 此 在 交叉 


必须 对 种 群 中 的 所 有 维 
产生 的 子 代 公式 如 下 : 


进行 归 


化 处 


里 。 每 次 纵向 交叉 操作 


不 
前 
后 


MS (i,d)=r:X(i,d)+(—r):X(i,d,) (6) 


i e N(1,M ),d,,d; e N(,D),r el0,1] 


其 中 : MS,.(i,d1) 是 粒子 XQ 的 第 di 维和 第 d, 维 

得 到 的 子 代 。 纵 向 交叉 概率 p,，( 一 般 取 0.2 一 0.8) 
能 使 粒子 摆脱 局 部 最 优 的 状况 ， 交 叉 后 进入 竞争 机 制 ， 将 产生 
的 子 代 与 其 父 代 粒 子 比 较 ， 选 择优 秀 粒子 保留 在 DS,. 中 。 


进行 纵向 交 


交叉 操作 


4 ”纵横 交叉 算法 优化 神经 网 络 


BP 神经 网 络 09 是 指 基 于 产生 误差 反 向 传播 算法 的 前 馈 


经 网 络 , 它 是 由 输入 层 、 
隐 含 层 的 功能 是 连接 输入 和 输出 。 


| 


力 和 噪声 赛 限 ， 已 广泛 应 用 于 不 同 的 j 


网 络 在 预测 过 程 中 ， 采 | 


梯度 下 降 法 调整 权 值 和 闷 值 ， 经 党 


后 得 到 的 子 代称 为 中 庸 解 ( MS,。 ，MS,。 ) 。 所 得 到 的 子 代 与 其 
父 代 相 比 ， 保 留 适应 度 更 好 的 解 ， 称 为 占 优 解 〈 横 向 交叉 占 优 
解 DS,. ， 纵 向 交叉 占 优 解 DS,. )。 


.3.1 横向 交叉 操作 
> 横向 交叉 与 遗传 算法 中 的 交叉 相似 ， 横 向 交叉 (横向 交叉 
>< 概率 通常 取 1) 是 在 两 个 粒子 中 做 算数 交叉 操作 ， 且 两 粒子 是 
= 在 同一 维 随机 产生 。 与 传统 交叉 相 比较 ， 加 强 了 局 部 最 优 粒子 
三。 的 边缘 搜索 能 力 ， 增 大 了 搜索 的 范围 。 交 又 后 产生 的 子 代 公式 
OO 如 下 : 
MS (i,d)=nxX(i,d)+(—n)xX(,d)+ 
cx(X(i,d) -xX(j,d) 
MS (1,d)=r, xX(j,d)+(—r,)xX(i,d)+ 
c, x(X(],d)—X(i,d)) 
i, jeN(I,M), deN(d,D) 
其 中 : a、 是 [0 和 ] 的 随机 数 ，ci 、c; 是 [1 中 的 随机 数 ，M 


为 粒子 数 的 范围 ，D 为 变量 的 维 数 ， 关 (i,d) 、X(j,d) 各 自 代 
表 父 代 粒 子 站 (i 和 (站 的 第 4d 维 ; MSi.(i,d)、MSi.(jsd) 各 自 
代表 (i,4) 和 站 (j,d) 进行 横向 交叉 后 得 到 的 第 d 维 子 代 。 

以 社会 学 的 角度 来 看 ， 式 (4) (5) 的 第 一 项 是 粒子 此 时 最 好 
的 值 ， 第 二 项 表示 父 代 粒子 相互 作用 。 前 二 项 通过 nn (惯性 权 
重 因子 ) 较 好 的 结合 在 一 起 。 式 (4)(5) 中 的 第 三 项 以 较 小 的 概率 
在 种 群 边缘 寻 优 ， 增 大 了 搜索 范围 。 得 到 的 子 代 将 与 父 代 相 比 
较 ， 选 择 适应 度 更 好 的 保留 在 占 优 解 DS;. 中 ， 
代 。 


进行 下 一 次 达 


的 预测 能 力 ， 如 图 


线 , 有 昌 CSO-BP、PSO-BP 两 种 模型 


一 致 ， 实 验 数据 都 是 来 电力 市 场 1 


出 现 陷入 局 部 最 优 的 问题 ， 随 着 干 
时 ， 神 经 网 络 参 数 计算 量 也 会 急剧 
而 CSO 算法 全 局 搜索 效果 较 好 ， 能 很 大 程度 上 提高 神经 网 
7 所 示 ， 是 两 种 模型 的 1 步 预测 最 优 收敛 


扰 因 素 和 样本 数据 逐渐 增 
增加 ， 收 敛 速度 也 将 变 缓 。 


又 
后 


神 


名 含 层 、 输 出 层 组 成 的 三 层 网 络 结构 ， 


由 于 其 具有 良好 的 自学 习 能 
预测 领域 。 由 于 BP 神经 


人 
EA 


加 


络 


的 种 群 大 小 和 和 迭 代 次 数 设 


o 


假设 


m ,变量 维 数 D=nxh 


mxh+h+m 


经 网 络 输入 层 , 隐 含 层 和 输出 层 的 节点 数 是 a 、h 


置 


、 


如 下 均 方 误差 公式 (适应 度 函 数 ) 来 获得 : 


N m 


f= (ppb,) 


w=1 1=1 


且 


输出 ，NN 表示 训练 样本 数 。 
CSO 优化 凶 


上 ， 六 表示 神经 网 络 的 


中 :jp, 表 示 神 经 网 络 的 实际 输 昌 


入 


。 每 个 粒子 的 适应 值 使 用 


经 网 络 的 过 程 概括 如 下 : 


a) 按 照 本 文 实验 的 训练 样本 ， 确 定 BP 的 层 数 和 各 层 的 


经 元 数目 ， 并 进行 实数 编码 。 


林 、 


神 


b) 确 定 本 文中 主要 参数 :种 群 规模 、 最 大 迭代 次 数 、 纵 向 交 
又 概率 p, 。 在 编码 的 解 空 间 中 ， 随 机 生成 初始 种 群 碟 。 

9) 结 合式 (7) 计 算 群 体 中 所 有 粒子 的 适应 值 。 

dd) 按照 式 (4)(5) 进 行 横向 交叉 操作 , 获得 的 子 代 保存 在 矩阵 


MS 里 


看 ， 再 计算 该 矩阵 里 所 
与 父 代 素 《〈 即 DS,.， 第 一 代 除 儿 
子 保留 在 DS. 中。 


e) 根 据 式 (6) 进 行 纵向 交叉 后 7 
然后 计算 MS,。 每 个 粒子 的 适应 值 ， 
选择 更 优秀 的 粒子 保留 在 DS,. 中 。 


生 的 解 保存 在 矩阵 MS, 


有 粒子 的 适应 值 , 将 得 到 的 结果 
) 对 比 ， 选 出 适应 度 更 好 的 粒 


与 其 父 代 头 ( 即 DS,. ) 对 上 


kk 
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了 判断 是 否 达到 设 定 的 要 求 。 如 果 和 迭代 次 数 大 于 所 设 定 的 
最 大 迭代 值 ， 则 和 迭代 结束 。 并 将 DSy 中 适应 度 最 好 的 一 组 解 设 
为 为 神经 网 络 所 对 应 的 权 值 和 阔 值 。 否 则 ， 转 到 步骤 四 再 一 次 
友 代 。 


10 20 30 40 50 60 70 80 90 100 


迭代 次 数 


图 7 不 同 算法 最 优 收敛 曲线 


5 ” 算 例 分 析 


为 了 验证 本 文 所 提 模 型 EEMD-VMD-CSO-BP 在 一 步 和 多 
步 电 价 预 测 中 的 有 效 性 ， 因 此 进行 了 以 下 二 项 实验 研究 加 以 证 


解 技术 最 好 , 为 了 进一步 分 析 本 文 提出 的 二 层 分 解 技术 和 CSO 
优化 算法 的 混合 模型 所 具有 的 优势 ,在 本 项 实验 中 进行 了 以 下 
三 类 比较 ， 对 比 了 不 同 模型 在 单 步 预 测 中 误差 结果 ， 如 表 2 所 
示 。 


ne 


一 类 比较 : VMD-CSO-BP 模型 和 EEMD -CSO-BP 模型 
相 比 CSO-BP 模型 平均 百分比 误差 分 别 减少 了 46% 和 62%, 结 
果 显 示 了 单 分 解 技 术 的 优势 。 

第 二 类 比较 : EEMD-VMD-CSO-BP 模型 相 比 VMD-CSO- 
BP 模型 和 EEMD-CSO-BP 模型 平均 百分比 误差 分 别 减少 了 51% 
和 30%, 所 以 本 文 提出 的 二 层 分 解 技 术 优 于 单 分 解 技术 , 此 外 ， 
图 8-10 可 以 看 出 , 随 着 步 长 的 增加 , 二 层 分 解 技 术 对 预测 精 
度 影响 越 来 越 小 ， 原 因 是 由 于 步 长 的 增加 ， 除 IMF1 外 的 其 他 
IMEF 形成 了 误差 的 累积 ， 因 此 在 第 二 层 分 解 中 VMD 对 提高 预 
测 精度 的 作用 会 越 来 越 大 。 
第 三 类 比较 :CSO-BP 模型 较 BP-NN 模型 的 平均 


人 


百分比 误 
差 减 少 了 11%: 说 明 CSO 算法 对 神经 网 络 产生 了 积极 的 影响 ， 

通过 优化 BP 的 参数 ， 使 其 产生 接近 最 优 的 权 值 和 阔 值 ， 从 而 
有 更 高 的 预测 精度 。 
在 第 二 项 实验 研究 中 ， 所 有 预测 结果 如 表 3 所 示 ， 显 然 本 
文 模型 在 一 步 三 步 五 步 中 预测 误差 最 小 ， 且 分 别 在 不 同 的 电力 
市 场 都 适用 ,进一步 证 明了 所 提 模 型 在 预测 能 力 上 较 其 他 混合 


明 ,， 所 有 模型 的 迭代 次 数 设 为 100, 种 群 大 小 M 设 为 20， 且 神 
经 网 络 参数 均 设置 相同 。 在 第 一 项 实验 中 ， 以 新 南 威 尔 士 州 的 
电价 数据 作为 实验 数据 ,使 用 所 提 模 型 的 组 成 部 分 : BP、CSO- 
BP、VMD-CSO-BP 和 EEMD-CSO-BP 被 选 为 基准 模型 ， 与 本 
文 模型 进行 对 比 ， 验 证 了 混合 策略 的 性 能 ， 在 第 二 项 实验 中 ， 
为 了 全 面 评估 所 提 模 型 ， 本 文 分 别 以 新 南 威尔士 州 和 昆 士 南 州 
电价 为 实验 数据 , 将 三 个 现 有 模型 , 包括 WT-GA-SVM、EMD- 
GA-BP 和 PSO-BP 与 本 文 模型 进行 比较 。 另 外 ,， 采 用 了 三 种 标 
准 来 评估 所 有 模型 的 预测 效果 ， 如 式 8、 式 9、 式 10 分 别 是 平 
均 绝 对 误差 (mean absolute error，MAE ) 、 均 方 根 误差 (root 
mean square error，RMSE) 以 及 MAPE 的 表达 式 。 


Ma 二 20- 交 (8) 


RMSE = 0-x0) (9) 


1=1 


_ Lx(O-x(). 

MAPE We x 0 (10) 

其 中 ; N 表示 测试 的 数据 数目 ，X(?) 是 指 在 时 间 为 t 时 候 的 实 
际 电价 值 ， 豫 () 是 指 在 t 时 刻 的 预测 电价 值 。 


表 1 显 示 了 本 文 模型 和 基准 模型 在 一 步 和 多 步 预测 中 的 预 
测 误差 , 显然 本 文 模 型 的 MAE、RMSE、MAPE 在 一 步 (S1) 、 
三 步 (S2)、 五 步 (S3) 中 的 值 最 小 , 所 示 结 果 证 明了 本 文 模型 的 分 


模型 具有 明显 的 优势 。 
表 1 不 同 模型 的 预测 性 能 分 析 电 力 市 场 1) 
预测 步 长 1 步 3 步 5 步 
BP-NN 预测 模型 
MAE 24.54 31.23 31.45 
RMSE 40.99 46.23 47.77 
MAPE(%) 29.43 45.67 47.86 


CSO-BP-NN 预测 模型 


MAE 22.51 29.81 30.01 

RMSE 37.75 40.34 42.66 

MAPE(%) 26.85 41.76 44.42 
EEMD-CSO-BP 预测 模型 

MAE 8.62 12.82 15.25 

RMSE 15.58 22.45 26.56 

MAPE(%) 10.11 18.90 24.21 
VMD-CSO-BP 预测 模型 

MAE 7.81 8.22 9.05 

RMSE 10.69 11.99 13.30 

MAPE(%) 14.58 15.01 15.88 

EEMD-VMD-CSO-BP 预测 模型 

MAE 3.97 5.55 7.69 

RMSE 5.60 7.76 10.34 

MAPE(%) 7.04 9.62 14.85 
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表 2 表 1 中 的 单 步 预测 结果 对 比 (电力 市 场 1 结果 ， 可 以 发 现 二 层 分 解 技术 在 提高 模型 的 预测 能 力 上 优 于 单 
模型 MAE(%) RMSE(%) MAPE(%) 分 解 模 型 ， CSO 算法 对 BP 模型 产生 了 积极 影响 ， 通过 优化 其 
nD CB 输入 层 和 隐 含 层 的 权 值 和 闵 值 ， 从 而 提高 预测 精度 。 与 本 文 所 
a 49% 53% 5s1% 提 其 他 混合 模型 相 比 ，EEMD-VMD-CSO-BP 混合 模型 具有 最 
MD CUO BE 佳 性 能 ， 表 明 所 提 模 型 非常 适用 于 非 稳定 性 电价 预测 。 
EEMD-VMD-CSO-BP 
300 
VS. 54% 64% 30% 一 一 实际 值 
一 一 BP-NN 预 测 值 
EEMD-CSO-BP 250 上 CS0-BF- 预测 信 
WU，，| 一 和 一 EEMD-CS0-NN 预 测 值 
VMD-CSO-BP 插 | 一 "一 VMD-CS0-NN 预 测 值 
200 一 一 EEMD-VWMD-CSO-NN 预 测 值 
VS. 65% 72% 46% s 
z 
CSO-BP-NN 金 150 
J 
EEMD-CSO-BP 型 前 
VS. 62% 58% 62% 
CSO-BP-NN 0 
CSO-BP-NN 0 1 1 1 | 1 1 1 1 1 
600 610 620 630 640 650 660 670 680 690 700 
VS. 8% 10% 11% 时 间 ( 半 小 时 ) 
| 图 8 提前 一 步 电价 预测 结果 (电力 市 场 1) 
3 ”本 文 模型 与 已 存在 模型 的 预测 性 能 对 比 -i i 
一 一 实际 值 
EEMD-VMD-CSO-BP WT-GA-SVM 一 = 一 BP-NN 预 测 值 
一 一 C80-BP-NN 预 测 值 


一 一 EEMD-CSO-NN 预 测 值 
一 一 VD-CS0-NN 预 测 值 
一 + 一 EEMD-VMD-CS0-NN 预 测 什 


Sl S3 S5 Sl S3 S5 


MAE 3.97 333 7.69 10.28 15.67 20.63 200 


RMSE 5.60 7.76 10.34 16.42 22.49 33.38 


MAPE 7.04 9.62 15.01 13.25 19.31 26.63 


电价 ($/MWh) 
图 


己 
S 


电力 市 场 EMD-GA-BP PSO-BP 


1 Sl S3 S5 Sl S3 S5 EL) y 


MAE 18.56 25.79 30.21 21.17 34.26 45.96 


0 
600 610 620 630 640 650 660 670 680 690 700 


RMSE 20.59 33.77 38.25 30.52 42.52 55.95 时 间 ( 半 小 时 ) 
MAPE 15.50 21.92 29.96 24.15 32.37 41.69 图 9 ”提前 三 步 电 价 预测 结果 (电力 市 场 1) 


EEMD-VMD-CSO-BP WT-GA-SVM 
Sl S3 S5 Sl S3 S5 


一 一 实际 人 
MAE 3.27 5.29 8.68 9.51 13.29 19.53 加 
250 上 一 一 CS0-BP-MN 预 测 值 
一 一 EEMD-CSO-NN 预 测 值 
一 + 一 VWD-CS0-NN 预 测 什 
一 # 一 BEMD-VMD-CS0-NN 预 测 值 


RMSE 4.73 7.48 12.05 13.93 20.41 30.64 


MAPE 5.48 8.50 13.88 11.22 18.67 24.21 会 200 

电力 市 场 EMD-GA-BP PSO-BP 2 
她 
2 Sl S3 S5 Sl S3 S5 


己 
[=4 


MAE 14.44 21.56 27.24 19.78 33.86 42.35 
RMSE 16.51 28.78 34.76 28.60 40.55 53.44 50 nak / 


MAPE 12.20 18.37 27.13 22.99 30.86 40.57 


6 ”结束 语 图 10 ”提前 五 步 电 价 预测 结果 (电力 市 场 1) 
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